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Abstract 

By explicit construction, we show that one can in a simple way introduce and 
measure gravitational holonomies and Wilson loops in lattice formulations of 
nonperturbative quantum gravity based on (Causal) Dynamical Triangulations. 
We use this set-up to investigate a class of Wilson line observables associated with 
the world line of a point particle coupled to quantum gravity, and deduce from 
their expectation values that the underlying holonomies cover the group manifold 
of SO(4) uniformly. 
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1 Introduction 


Our quest for a theory of quantum gravity meets with numerous challenges. We 
not only have to define the theory nonperturbatively, but must also ascertain that 
it actually exists and has desirable physical properties, including a well-dehned 
classical limit. Many steps in this construction must be formulated in terms of 
observables, which in a diffeomorphism-invariant theory are notoriously hard to 
come by. In this article, we focus on a particular class of observables, involving 
gravitational Wilson loops, and a particular candidate theory of nonperturba- 
tive quantum gravity. Causal Dynamical Triangulations (CDT) [1, 2, 3]. In this 
theory, the challenges mentioned above pose themselves in very concrete terms 
and can also be addressed concretely, including the use of powerful numerical 
methods. 

To understand our analysis of Wilson loops, a comprehensive understand¬ 
ing of CDT quantum gravity will not be necessary^; we will conhne ourselves 
to a brief description of the approach, and in later sections give some details of 
the geometric set-up, to the extent they are needed. In a nutshell, CDT is a 
covariant, quantum held-theoretic lattice formulation of gravity, where the non- 
perturbative sum over spacetime geometries is realized in terms of piecewise flat 
four-geometries. They are assembled from four-dimensional Lorentzian build¬ 
ing blocks in such a way that only causally well-behaved spacetime histories are 
included in the path integral. 

To perform the actual sum over these histories one must rotate them to Eu¬ 
clidean signature. It is important to understand that not all Euclidean triangu¬ 
lations lie in the image of the Wick rotation map, but only those in a subset, 
which carry a memory of the causal properties of their Lorentzian origin, most 
notably, the absence of “baby universes” (topology changes of spatial slices). The 
triangular building blocks or four-simplices are characterized by their side length 
a, which plays the role of a UV cut-off. The continuum limit of the regularized 
path integral involves a limit a —)■ 0, possibly accompanied by a readjustment 
of the bare coupling constants, such that physics stays invariant. In recent work 
[6] we demonstrated explicitly how a renormalization group flow is implemented 
in CDT quantum gravity, despite the absence of a background metric and the 
absence of any obvious correlation length. Apart from being a rather remarkable 
result, our analysis highlighted the need for further observables to provide inde¬ 
pendent checks on our condition of keeping physics constant while altering the 
renormalization group scale. 

Only a few observables are known in CDT quantum gravity and have been 
investigated quantitatively, including the volume prohle of the dynamically gen¬ 
erated quantum universe [3, 7], as well as its Hausdorff and spectral dimensions 

^The interested reader may consult our overview and review articles [4, 5]. 
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[2, 8 ]. Note that all of them involve measurements of lengths and volumes. This 
is in contrast with the classical continuum theory, where one describes the non¬ 
trivial, local structure of spacetime in terms of its curvature, which is a function 
of the derivatives of the spacetime metric g^u{,x). A key question we would like 
to answer is whether there is a meaningful notion of “curvature” or “quantum 
curvature” in nonperturbative quantum gravity, which on the Planck scale is well- 
dehned and yields hnite values, and on macroscopic scales goes over into one of 
the standard curvatures of general relativity. 

In Regge Calculus [9] and Dynamical Triangulations^ there is a simple, dis¬ 
cretely dehned expression for the local scalar curvature in terms of dehcit angles, 
which we will review in Sec. 2 below, but unfortunately it becomes singular in a 
naive continuum limit. This is not at all surprising since the continuum dehni- 
tion of the curvature involves second derivatives of the metric and a typical held 
conhguration g^u{,x) in the path integral is not expected to have well-behaved 
derivatives. 

A main motivation for considering Wilson loops to try to dehne some coarse¬ 
grained measure of curvature comes from gauge held theory. Here one can con¬ 
struct a nonlocal, gauge-invariant observable by taking the (trace of the) path- 
ordered exponential of the gauge potential A^{x) along a closed curve 7 , to obtain 
the so-called Wilson loop [11] 

W^{A) =TrVexp ([ A, (1) 

J7 

with V denoting path ordering. The relation with the local curvature tensor 
F^u{x) is exhibited by expanding the path-ordered exponential (the holonomy) 
around an inhnitesimal square loop of side length e in the /iz/-plane, yielding 

Vexp^A = l+gF;;,Xae^ + 0{e^), ( 2 ) 

where Xa are the generators of the Lie algebra of the gauge group and g denotes 
the coupling constant. Moreover, the scaling behaviour of large Wilson loops 
provides a test for whether the theory is conhning. Wilson loop observables are 
robust in the sense that they have a natural representation in terms of lattice 
variables in lattice gauge theory and have been used successfully in numerical 
studies. 

In gravity, one can use the metric-compatible Levi-Civita connection r^^{x) 
to construct holonomies and gravitational analogues of Wilson loops, as we will 

^Lattice gravity in terms of Dynamical Triangulations (DT) is the purely Euclidean precursor 
of Causal Dynamical Triangulations, see, for example, [10]. Our theoretical considerations about 
Wilson loops presented below, up to and including Sec. 5, are also applicable to DT. 
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describe in more detail in Sec. 2 below. The path-ordered exponential of T along 
a path dehnes a notion of parallel transport of tangent vectors, and all physical 
information contained in the Riemann curvature tensor Rx^ui^) can be retrieved 
from suitable inhnitesimal holonomies, analogous to the situation in gauge theory 
captured by eq. (2). However, Wilson loops are not diffeomorphism-invariant, 
unless the underlying loops are dehned in physical terms. Of course, this does not 
mean that one cannot construct quantum observables that depend on holonomies 
or Wilson loops and are diffeomorphism-invariant. 

Gravitational Wilson loops on spacetime have been little studied, with the 
exception of work in perturbative quantum gravity [12] and in the context of the 
search for a nonabelian Stokes’ theorem [13]. The story is different in canonical 
quantum gravity, where holonomies along spatial curves play a prominent role in 
Loop Quantum Gravity [14]. This approach differs radically from perturbative 
quantum gravity where the dynamical variables are local helds like the metric 
g^iy{x). Instead, in loop quantum gravity nonlocal holonomies are taken as part 
of a set of fundamental variables in terms of which the entire quantum dynamics 
should be expressed. In the quantum theory they are promoted to hnite opera¬ 
tors, which are assumed to not need any renormalization. This is different from 
ordinary gauge theory, where the expectation values of Wilson loops need to be 
renormalized. 

In this article we consider quantum gravity in the GDT formulation. Despite 
being nonperturbative, it is nevertheless an ordinary quantum field-theoretical 
framework. In order to extract physical information from suitable loop averages 
when the lattice cut-off is taken to zero, we therefore expect that observables 
involving Wilson loops will require renormalization. 

Motivated by the fact that Wilson loops - at least inhnitesimal ones - encode 
retrievable curvature information, and encouraged by their success as observables 
in nonperturbative QGD, our ultimate goal is to construct and measure quantum 
curvature observables in nonperturbative quantum gravity based on holonomies 
or Wilson loops. As explained earlier, they should also provide us with a no¬ 
tion of averaging or coarse-graining^, to allow for a comparison with ordinary 
macroscopic curvature in a semiclassical limit. We do not know a priori whether 
such observables exist, and we are not aware of an explicit construction in any 
approach to nonperturbative quantum gravity.^ The results derived in this pa- 

•^We note in passing that the analogous averaging problem in classical general relativity has 
not been resolved (see, for example, [15]). 

■^Theoretical arguments were put forward in a different formulation of lattice gravity based 
on Regge calculus, promoting large Wilson loops as carriers of nonperturbative information 
[16]. Although sympathetic to the aim, we are unable to follow the technical claims in [16] 
or to understand how the construction can be implemented meaningfully in a nonperturbative 
context. 
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per, involving both theoretical considerations and numerical simulations in four 
dimensions, hopefully present a step in the direction of our main goal, as well as 
demonstrating that CDT quantum gravity as a framework is perfectly suited to 
studying observables of Wilson loop type. 

In what follows, we begin by reviewing holonomies in continuum gravity (Sec. 
2 ), as well as their counterparts in piecewise flat spaces and, more specihcally, 
in dynamical triangulations (Sec. 3). In Sec. 4, we introduce the invariant an¬ 
gles characterizing a general SO(4)-holonomy, and derive an explicit expression 
for an associated distribution of their possible values on the group manifold. A 
convenient choice of coordinate frames on the four-simplices of the triangulations 
is introduced in Sec. 5, as well as two different ways to compute the holonomies 
of closed lattice loops. In Sec. 6 the discussion focusses on a specihc class of Wil¬ 
son loops, associated with the world line of a point particle, and their concrete 
implementation in the full, nonperturbative CDT path integral. After a brief de¬ 
scription of the Monte Carlo simulation of the combined gravity-particle system. 
Sec. 7 contains our main computational result, the measured distribution of the 
invariant angles for the class of Wilson lines considered. We conclude in Sec. 8 
with a discussion and outlook. 


2 Holonomies in gravity 


The Levi-Civita connection r((,^(a;) of a Riemannian manifold M with metric 
dfiuix) dehnes a notion of parallel transport of a vector along a curve 7 ^(A). 
Transporting along the curve between parameter values Ai and Af results 
in a general linear transformation of the vector, which is given in terms of the 
path-ordered integral, 

= ,V‘-(xd, (r.,)7 = r"„, (3) 


where V denotes path-ordering, the dot indicates differentiation with respect to 
the path parameter A, and = 7 (Ai) and X{ = 7 (Af) are the initial and hnal 
points of the path in M. 

Under a coordinate transformation x —)■ x{x), with M^^{x) = the 

path-ordered integral transforms nontrivially at its two endpoints Xi and Xf, 









( 4 ) 


in accordance with the transformation behaviour of vectors under coordinate 
transformations, namely, 

U^(x) = M^^(a;) W(a;). (5) 
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Before turning to the case of a piecewise flat manifold, let us look at the concrete 
construction of the path-ordered product Pexp(—J^F) for a given curve 7 (A). 
In general, 7 (A) will pass through several coordinate patches t/fc, k = 1,... ,n, 
with corresponding coordinates x^. 

Let us consider the simplest situation of this kind, where the initial point x\ 
of the curve lies in an open neighbourhood Uq and the hnal point Xf in an open 
neighbourhood 17i, such that Xi ^ t/i, Xf ^ Uq and the intersection Uq fl Ui is 
not empty. To perform the path integration along 7 of the connection F, an 
intermediate point 7(Amid) = a^mid must be chosen in the overlap region Uq fl f/i, 
and the integration performed in two pieces: from Ai to Amid over the connection 
Fo(xo) in the coordinates Xq of patch f/g, and subsequently over the connection 
Fi(x5‘) in terms of the coordinates x^ of Ui. In addition, to account for the 
change of coordinate system, a matrix M(xmid) = to be inserted at 

the midpoint, leading to a combined expression schematically given by 

(Pe- M(xmid)"A(Pe- (6) 

Using the transformation law (4), it is straightforward to show that the value of 
expression ( 6 ) is independent of the choice of midpoint Xmid ^ UqP\Ui. 



Figure 1: Starting at an initial point xi, a path 7 (A) traverses a sequence of 
coordinate neighbourhoods f/fc, ending up at Xf. A path-ordered integral along 7 
of a connection is computed piecewise in every patch [Uk, The switch from 

the to the {k + 1 )®* coordinate system can happen at an arbitrary midpoint 
Xmid,k ^ Ukr\ Uk +1 and is associated with a matrix M(xmid,A:) as described in the 
text, c.f. eq. ( 8 ). 

If the path 7 runs through several coordinate neighbourhoods, the construc¬ 
tion for the two patches just given can be reiterated, see Fig. 1. We are specihcally 
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interested in closed paths and therefore will consider the situation where 7 starts 
at Xi in neighbourhood (I/ojIxq}), passes through a sequence of n neighbour¬ 
hoods {Uki {Xk})) A; = 1 ,..., n, via their non-empty intersections Uk H Uk-i 7 ^ 0 , 
and hnally into Un+i = Uo, where it ends up at the same point X{ = Xi it started 
from. The path-ordered integral or holonomy associated with the oriented loop 
7 based at Xi is then represented by 

n 

fc=l 

where it is understood that the matrix multiplication is from right to left as 
the loop parameter A in 7 (A) increases from Ai to Af. While the path-ordered 
integral appearing in eq. (3) is valued in GL(4,1R), the holonomy matrix (7) 
on an orientable manifold is valued in SO(4). The transformation matrix from 
coordinates {x^_{\ to {x^} is given by 

M{x^,,kT^ = ( 8 ) 

and is evaluated at the midpoint Xmid,?;, to be chosen freely along 7 in the 
overlap region Uk O Uk-i- The holonomy (7) still depends on the initial or base 
point Xi and under a coordinate transformation x —)■ x{x) will transform according 
to (4) as 

= M{xi) (Ve~ M~^{xi), with M^y{x) = (9) 

It follows that the conjugacy class of the holonomy matrix 'Pexp(—^^T) is co¬ 
ordinate independent. It is also easy to show that it does not depend on the 
starting point Xi chosen along the loop 7 (A). In this paper we will precisely study 
such ccordinate-independent conjugacy classes. 

3 Holonomies in Dynamical Triangulations 

Studying holonomies in the context of piecewise flat geometries simplihes the 
above discussion considerably, as we will see. The building blocks of the piece- 
wise linear geometries used in dynamical triangulations are identical, equilateral 

®In causal dynamical triangulations one usually works with two different edge lengths, one 
for time-like and one for space-like edges [5]. In the present study we consider for simplicity 
the special case where after the Wick rotation of CDT the two edge lengths are identical, and 
each triangulation therefore becomes equilateral. 
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four-simplices, which by assumption are everywhere flat on the inside, like the 
building blocks of Regge calculus [ 17 ]. When these building blocks are glued 
together along identical boundary three-simplices or “faces” to construct a four¬ 
dimensional piecewise flat manifold, curvature will generically appear in a singular 
fashion along the two-dimensional subsimplices of the triangulation, the triangles 
or “hinges”. 

Recall that the geometry of a four-simplex is completely fixed by its edge 
lengths (in our case the single edge length £), and that the geometric proper¬ 
ties of a four-geometry assembled from such simplices are encoded in the gluing 
data (how faces are identified pairwise), neither of which requires the introduc¬ 
tion of coordinates. Indeed, an important strength of the nonperturbative path 
integral formulation of CDT comes from the fact that no coordinates have to be 
introduced, and that the path integration does not contain unphysical coordinate 
reparameterizations or other parameter redundancies. 

In the present piece of work we are not going to change the way we perform 
the path integral, but in order to analyze particular quantum operators involving 
holonomies we will introduce coordinate systems on individual four-simplices. 
We are in principle completely free how to do this. Since the final result will not 
depend on these choices, it is convenient to use the same Euclidean flat coordinate 
system on every simplex. (We will specify our particular choice later on.) Since 
this makes the metric constant, the connection E vanishes everywhere on the 
four-simplex, and its path-ordered integral along any curve 7 is the unit matrix, 
as long as 7 remains inside the simplex. 

Our considerations about computing holonomies of loops passing through sev¬ 
eral coordinate patches apply to these simplicial manifolds as follows. Two neigh¬ 
bouring four-simplices si and S2 with associated flat Euclidean coordinate systems 
{xi} and {X2} always have a three-simplex (tetrahedron) a in common. Since 
there is no curvature associated with (the interior of) this three-dimensional in¬ 
terface, combining it with the interiors of the two four-simplices results in a single 
open set whose geometry is flat and constant. Calling Ui the coordinate patch 
parameterising (the interior of) four-simplex Si, this implies that (Ri,{a^i}) can 
be continued to the interior of S2 and/or {U2, {x^}) can be continued to the in¬ 
terior of Si to create a non-empty overlap region Ui fl U2- Like in the smooth 
case above, this allows us to associate a path-ordered integral with any curve 
passing from si to S2 by inserting a matrix M = associated with the change 
of coordinates xi —)■ X2{xi) in between the piece of the path-ordered integral 
computed in Ui and that computed in U2. Since the coordinate systems are flat 
Euclidean and moreover are the same for all simplices, the transition matrix M 
on the four-dimensional overlap region Ui (11/2 is constant and given by a four¬ 
dimensional rotation. Reverting once again to a simplicial description, we may 
therefore simply associate the matrix M with the entire (interior of the) three- 



dimensional interface a between si and S2. To capture this simple dependence, 
we will introduce a new notation for the corresponding rotation matrix, namely, 

=: R{s 2 , SiTu e 0 ( 4 ), ( 10 ) 

with the implicit understanding that R still depends on the coordinates {x^} 
and {X2}. The explicit form of R{s2,Si) depends on the relative orientation of 
the two coordinate frames (si, {x^}) and (s2, {2^2})) be computed once 

the coordinate systems have been specihed. We conclude that the path-ordered 
integral along any path crossing from Si to S2 anywhere in the interior of a will 
pick up a factor of R{s2, si). 

Consider now a closed path 7(A) in a piecewise flat simplicial manifold T, and 
assume that it does not pass through any of the two-dimensional subsimplices 
of T, thereby avoiding potential curvature singularities. The holonomy along 
7 is then given by the ordered product of the rotation matrices associated with 
subsequent crossings of 7 from one four-simplex to the next. Since the connection 
inside the four-simplices vanishes, we can restrict ourselves to a limited set of 
standardized closed paths without losing any holonomy information. For our 
present purposes it is convenient to use only loops consisting of straight segments 
between the centres of neighbouring four-simplices. For the holonomy Rl of such 
a loop L, which passes through a sequence si, S2, ■ ■ ■, s„, si of four-simplices, 
formula ( 7 ) reduces to a product of the corresponding rotation matrices, 

Rl 7?(si, R(^Snj >s^_i) ■ ■ ■ 7?(s2, <si). (H) 

Like the general holonomy ( 7 ), Rl still has a residual coordinate dependence 
and transforms non-trivially at its base point under a coordinate transformation 
X —)■ x{x) on Si. Since we have already hxed the coordinate frames of four- 
simplices to be flat and Euclidean, such a coordinate transformation must be a 
0 ( 4 )-rotation A, and Rl will transform by conjugation accordingly, 

R, AR,AL ( 12 ) 

If we parallel-transport a vector V in the tangent space to (si, {xi}) around the 
loop L, it will undergo a four-dimensional rotation to a new tangent vector 


El = RlV. 


The angle between the original and the rotated vector, dehned as 


^v,Vl arccos 




( 13 ) 


( 14 ) 
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Figure 2: Local curvature associated with a hinge r of codimension two in a d- 
dimensional simplicial manifold: a vertex in d = 2, an edge in (i = 3 and a triangle 
in (i = 4. Parallel transport of a vector perpendicular to a hinge r along a minimal 
closed path -L('r) through the centres of the d-simplices which form the star of r 
rotates the vector by an angle equal to the dehcit angle associated with r. (For 
ease of presentation, for d = A the four-dimensional star of r is not shown.) 


is independent of A, since a S'0(4)-rotation preserves scalar products V ■ W oi 
vectors in IR^. 

The set-up we have just introduced is closely related to how local curvature 
is described in Regge calculus. Generalizing from four dimensions, the curvature 
of a d-dimensional simplicial manifold is located at the subsimplices of dimension 
d — 2. Given such a hinge r, which does not lie on the boundary of the manifold, 
the curvature associated to it can be thought of as the Gaussian curvature of 
a small two-dimensional surface perpendicular to r. The surface is locally flat, 
with a conical singularity at the location where the hinge meets the surface. The 
magnitude of the curvature is quantihed in terms of a dehcit angle, which in a 
given simplicial manifold can be extracted via parallel transport as we will explain 
below for the case d = 4. 

To better understand the geometry of the situation, consider hrst the local 
conhguration of d-simplices sharing a hinge r. They form the so-called star of r, 
which topologically speaking is a d-dimensional ball. They also form a circular 
neighbourhood around r, in the sense that we can construct a minimal closed 
piecewise straight path L{t) encircling r, which connects the centres of adjacent 
d-simplices in the star (see Fig. 2). 

Returning to the physically relevant case d = 4, how will an arbitrary vector V 
be affected by parallel transport along L{t)7 Note hrst that V will nothe rotated 
at all if in the initial four-simplex si it is parallel to the two-plane spanned by 
the triangle r (identifying the linear structure of the simplex and its tangent 
space). This can be easily understood as follows. In si, choose an orthonormal 
coordinate system {x^} such that the hinge r lies in the plane spanned by xf 
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and x\, say. Since r is common to all four-simplices in the star of r, one 
can in each of them make the same choice for the two coordinate axes and 
x\. relative to r. In other words, an arbitrary vector = (0,0, in si 

will have exactly the same form in each of the Sk, independent of the choice 
of the coordinates in the directions perpendicular to r. It follows that only 
components of V perpendicular to r can be affected nontrivially by the holonomy 
matrix. Parallel transport around r along the minimal loop L(t) will therefore 
map a vector = {V^,V‘^,0,0) to some = (V^, V^,0,0). The only S'0(4)- 
transformations that can have this effect belong to the one-parameter SO{2)- 
subgroup mapping the plane orthogonal to r into itself. As a consequence, the 
rotation undergone by V is characterized by a single angular parameter e, where 


cose = , — — , for 1/ T r. (15) 

W^^VLiryVLir) 

We recognize e as the dehcit angle associated with the triangle r, as dehned in 
Regge calculus, and note its coordinate-independent character. 

To understand the range of the coordinate-invariant information that can 
be obtained by studying holonomies, let us recall some properties of the group 
S'0(4). A maximal torus of S'0(4) is given by the two-parameter set of matrices 


17(^1, ^2) 


/ cos 6*1 
— sin 6*1 
0 

V 0 


sin 6*1 
cos 6*1 


0 

0 


0 0 \ 

0 0 

cos 6*2 sin 6*2 

— sin 6*2 cos 6*2 / 


( 16 ) 


forming a maximal abelian subgroup SO{2) x SO{2). This implies that the rank 
of the group is 2. Moreover, given a compact connected Lie group G ~ like 
- and a maximal torus H G G, each element g & G is conjugate to an element 
h E H, that is, there is a a; G G and a h E H such that g = xhx~^ (see, for 
example, [18]). 

By the action of the holonomy Rl{t) of a minimal loop L{t) around a triangle 
r, with the coordinate choice made above, a vector V orthogonal to the r-plane 
according to relation (15) will be rotated by an angle e, with corresponding 
holonomy matrix 


Rl{t) 


/ cose sine 0 0 \ 

— sin e cos e 0 0 
0 0 10 

\ 0 0 0 1 / 


( 17 ) 


Note that in the CDT setting e can only take one of a discrete set of values, 
because all building blocks are identical and therefore all length and angular 
variables describing them come in discrete units. 
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It is important to realize that the matrix (17) from the point of view of S'0(4) 
corresponds to a particular type of rotation, a so-called “simple rotation”, which 
for S'0(n)-rotations on IR”^ is dehned as a rotation that leaves a linear subspace 
of dimension n — 2 hxed. The rotation matrix (17), and any matrix obtained from 
it by conjugation, leaves a two-plane through the origin hxed, and is therefore 
an example of a simple rotation in S'0(4). In dimensions n = 2 and n = 3 every 
rotation is simple, but this is no longer the case for n = 4, where a generic rotation 
instead is characterized by two angles, and is conjugate to a matrix of the form 
(16), with both 6 *i 7 ^ 0 and 6*2 7 ^ 0. Geometrically, such a “double rotation” 
consists of two independent (and commuting) rotations in two two-planes which 
are mutually orthogonal and therefore share only one point, the origin, which is 
also the only point mapped into itself by this kind of rotation. 

The distinction between a simple and a double rotation for a 5'0(4)-holonomy 
matrix is intrinsic and independent of coordinates. There is no coordinate trans¬ 
formation which will convert a double rotation with 7 ^ 0 and 02 7 ^ 0 to a 
simple one. The fact that parallel transport around a single triangle r in a sim- 
plicial four-dimensional manifold T results in a simple rotation has to do with 
the nature of the curvature singularity located at r. Parallel transport around a 
more general loop in T will in general not lead to a simple rotation. The same 
holds for parallel transport around loops in a general curved continuum manifold, 
regardless of whether the loops are hnite or inhnitesimal.® 

4 Invariants from holonomies 

Rather than operating with equivalence classes of rotation matrices under con¬ 
jugation, a convenient way of extracting the coordinate-invariant information of 
a holonomy matrix Ri is to take its trace, Tr(i? 2 .)- Especially in the context of 
gauge held theory, where the path-ordered integral is taken over the local gauge 
connection, this quantity is known as a Wilson loop. Because of the cyclic prop¬ 
erty of the trace, it is invariant under conjugation, 

Tt{ARlA^)=Tt{Rl), AeSO{4), (18) 

which means that for a generic holonomy matrix R^ G S'0(4) we can from (16) 
dehne the invariant quantity 

tilRi) := ^ Tr(Ri) = cos 9i + cos 62 - (19) 

®The fact that two angles are necessary to characterize holonomy matrices up to conjugation, 
as soon as one considers non-minimal loops in a four-dimensional simplicial manifold T, seems 
to have been overlooked by the authors of [16]. This also holds when T is almost flat and 
holonomies do not deviate much from the identity matrix. 
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The fact that one can interchange the two (2 x 2)-blocks on the diagonal of the 
matrix (16) by an appropriate conjugation is reflected in the fact that the right- 
hand side of equation (19) is invariant under the exchange of 6 i and 62 - Assuming 
for the sake of dehniteness that the angles 6 i take values in the interval [0,27r], 
we can by conjugation achieve that ( 6 * 1 ,^ 2 ) t (27r — 6*i,27r — 62 ), which likewise 
leaves (19) invariant. 

To extract information about both angles 6 i separately, we can supplement 
expression (19) by a second invariant, 

^ 2 (-Rl) := ^ Tr(i?|) + 1 = cos^ 9i + cos^ 62 . (20) 

If we £x the range of the angles to [ 0 , 7 r[ and require 61 < 62 , say, the invariants 
ti and t 2 fix { 01 , 62 ) uniquely. 

In the four-dimensional simulations we will extract this coordinate-invariant 
information by measuring ti and t 2 for a variety of closed curves, and on vari¬ 
ous ensembles of simplicial CDT geometries. As pointed out earlier, the 0(4)- 
rotations which enter into the construction of these quantities are not arbitrary 
but belong to a discrete set of possible rotations between neighbouring simplices. 
We will investigate how this influences the measured invariants ti. More specifi¬ 
cally, we will extract from them the distribution of the angles Oi and compare it 
to the distribution one would obtain if the holonomy matrices were distributed 
uniformly over the group manifold of S'0(4). 

In order to do this, we need to derive the theoretical distribution of the Oi on 
S'0(4). Recall that these angles were introduced in the context of the maximal 
torus (16). They are two of a total of six parameters needed to label points of 
S'0(4). By the theorem quoted earlier, every g G S'0(4) can be obtained from 
an element of the maximal torus by conjugation. An explicit way of doing this, 
which introduces an explicit parametrization of the group manifold, is given by 


■= U(ipi,tp2)W(QJu2J2)U(ei,e2)W(QJi,U}2flHtpi,tp2)'^, 

( 21 ) 

with 0 < 9i < 2 n, 0 < Ui < n and 0 < (pi < 2 n, and where 


W{uJi,UJ2) = 


/ cos U 2 

0 

0 

sina;2\ 

0 

COSCUi 

sinwi 

0 

0 

1 

— sin oji 

coscai 

0 

sin ijj 2 

0 

0 

cos UJ 2 / 


( 22 ) 


A straightforward way to obtain the desired distribution P{di,d 2 ) of the Oi on 
S'0(4) is to compute the Haar measure in terms of {Xk} = {^ 1 , ^ 2 , 1 ^ 1 , 1 ^ 2 , 9 ^ 2 }, 

where we have adopted a collective notation Xki A; = 1 , 2 ,..., 6 for the six group 
parameters. We then integrate the associated volume form over the parameters 
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not in the maximal torus, resulting in a two-form p{di,d2)d6id62- Normalizing 
p{9i,92) then gives the distribution P(6*i,6*2). 

We obtain the Haar measure by hrst computing the left-invariant^ one-forms 
Lk = &lso known as Maurer-Cartan forms, which take values in the Lie 

algebra of S' 0 ( 4 ). A left- and right-invariant volume form on the group manifold 
is then given by 

a/ det Gki d^X: ( 23 ) 

which involves the “metric”® Qki constructed from the left-invariant one-forms 

= — Tr(LfcL/). ( 24 ) 


Explicitly, one hnds 


det Qki = 2® (cos 9 i — cos 6'2)^(cos(2a;i) — cos(2a;2))' 


= 2^^ sin^ 


9 i + 9 


2 \ . 4 

sm 


9i — 92 


sin^(a;i -|- UJ2) sin^(a;i — U2). ( 25 ) 


Since the determinant factorizes into a part depending on the torus variables 9 i 
and a rest, the same holds for the volume form ( 23 ) and we can simply read 
off p{9i,92) up to a proportionality constant. After normalization one hnds the 
searched-for distribution of the angles 9 i and 92 , 


m,^2) 



( 26 ) 


Below we will report on the measurement of the distribution P{ 9 i, ^2) by Monte 
Carlo simulations in CDT, and will hnd them to be in perfect agreement with 
formula ( 26 ). 


5 Implementation 

We will now show how Wilson loops can be computed in the set-up of (causal) 
dynamical triangulations. Despite the fact that we are in four dimensions, the 
process is entirely straightforward and easy to implement. As explained earlier, 
we work with CDT building blocks that after the Wick rotation are equilateral, 
and therefore are all identical with respect to their geometric properties. 

Of the many possible coordinate systems for an individual four-simplex s 
of this type, we choose an orthonormal frame whose origin coincides with the 

^Of course, using the right-invariant one-forms Rk = instead would lead to the same 

result. 

is symmetric and bilinear with non-negative eigenvalues, but has degeneracies at some 
points g G S'0(4). 
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barycentre of s, which is equidistant to the hve vertices of the simplex. For 
dehniteness, we fix a scale such that the vertices all have geodesic distance 1 
to the barycentre. This is a coordinate-independent statement, which implies 
that the edge length of the simplex is £ = a/5/2. After assigning labels Pi, 
i = 1,..., 5 to the hve vertices of s, we hx the coordinate system uniquely by 
choosing coordinates for the point Pi. Representing the four-tuple of coordinates 
of the point by a column vector T* = x{Pi), the explicit choice is 




/O 

0 

0 


VI 


0 0 7 —qX 


KJ 

0 

KJ 

P 

a 

I 

d 

2 

a 

' 

d 

2 

a 

a 

3 

3 

3 

1 

1 

1 

-iJ 

4 

4 

4 


with 



One easily verihes that the coordinate vectors satisfy 


x^ = l, ^Tj = 0 , Xi-Xj 
i 



7^ P 


(27) 


(28) 


(29) 


with respect to the usual scalar product of Cartesian coordinates on Euclidean 
space.® Our standard choice of coordinates on a simplex s will be dehned through 
relations (27), (28) and (29); given another simplex s' with vertices P', i = l,... ,5, 
and coordinates x', its vertex coordinates will therefore be the same. 


£'(P') = f(Pj). (30) 

Since the standard coordinates depend on a specihc labelling {P/ of the ver¬ 
tices, we are left with a residual, discrete coordinate freedom, associated with 
relabelling those vertices. The latter is given in terms of the permutation group 
5*5 of hve elements. Given such a permutation tt : i h-)■ there is an associated 

permutation of vertices Pj h-)■ P 7 r(j), which in tnrn corresponds to a linear trans¬ 
formation P^ of the coordinate system. Adhering to the column vector notation 
introduced earlier, is a (4 x 4)-matrix given by 

(^TTl ,^7T2 5 5 ^7V4 ) = Ptt • • (31) 

Since the four vectors Xi, i = 1 ,..., 4 are linearly independent, we can solve this 
equation to obtain 

P(^TTl 5 X.jY^ ) ■ {xi,X2,X3,X4) / (32) 

®This coordinate construction generalizes to d dimensions, where the d-l-1 vertex coordinate 
vectors of an equilateral d-simplex are required to satisfy the conditions (29), with the third 
relation substituted by A • aq = — / 
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Since all scalar products are preserved by virtue of the relations (29), this transfor¬ 
mation is necessarily orthogonal, V.,^ G 0(4)d° In other words, we have obtained 
a representation of the permutation group 5*5 in terms of orthogonal matrices in 
four dimensions. If the permutation is even, we have detVT, = l and Vt^ G S'0(4), 
if the permutation is odd, we have det'P 7 r = —1 and V.,, ^ S'0(4). 

Having made a coordinate choice for a given four-simplex and vertex labelling, 
we will now address the explicit construction and computation of the holonomy 
Rl associated with a closed loop L passing through a sequence of four-simplices, 
along an oriented, piecewise straight path through adjacent simplex (bary)centres, 
as described in Sec. 3 above. Note that a path that passes through a simplex s 
distinguishes a triangle in s, namely, the triangle shared by the two faces (tetra- 
hedra) through which the path enters and leaves s. One may view this triangle as 
the hinge around which the path “bends”. If the path follows a complete set of 
consecutive four-simplices which share a single interior triangle r of a triangula¬ 
tion, we are back to the situation of a minimal loop, whose associated holonomy 
~ up to conjugation - is a rotation matrix of the form (17). 

The task at hand is to construct the rotation matrices i?(sj+i, Sj) in the expres¬ 
sion for the holonomy ( 11 ) of a loop passing through the simplices Si, S 2 , ...,<§„, Si. 
Having hxed a standard coordinate system for a given labelling of the vertices, the 
remaining gauge freedom we have to compute Rl is how to pick the vertex labels 
for the members of the set {sj}. For illustration, we will consider two different 
ways of computing the rotation matrix R{si+i,Si) associated with two adjacent 
simplices Sj and Sj+i in the sequence. In the first one, the vertex labelling of sim¬ 
plex Sj+i is related to that of the previous simplex Si along the loop L, and in the 
second one, the labellings of Si and Sj+i are picked independently and arbitrarily 
beforehand. It is the latter we will use in the simulations later on. 

5.1 Choosing vertex labels (version 1) 

Let the vertices of s* be labeled by {1,2, 3,4,5}, with standard coordinates (27) 
assigned to them. Assume that the boundary tetrahedron shared by Si and Sj+i 
has vertex labels ( 1 , 2 ,3,4}, and that the distinguished triangle r around which 
the holonomy “rotates” has vertex labels {1,2,3}. Now choose the vertex labels 
in simplex Sj+i such that for the vertices shared by both four-simplices we have 

p. = p'.^ i = 1,2,3, P 4 = P 5 , (33) 

where the primed vertices refer to simplex Sj+i. Using our standard coordinates, 
this hxes uniquely a coordinate system {x'^} in Sj+i, and the vertex coordinates 

^°Recall that 0{n) has two connected components, SO{n) being the component connected 
to the identity. The determinant of any member of the other component is —1. 
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of the two simplices are related by 

f'(P,) =f(P,), * = 1,2,3, f'(P 4 ) =x(P 5 ). (34) 

(Fig. 3 illustrates the analogous situation in two dimensions.) To compute the 
matrix P(si+i, Si), we make use of the observation from Sec. 3 that the coordinate 
system {x^} thus dehned in Si naturally extends to the neighbouring simplex Sj+i, 
and vice versa for the coordinates The relation between the x- and the x'- 

coordinates of the same point P involves the rotation matrix R and a translation 
(c.f. Fig. 3) according to 


S (P) = Si) x(P) + £, 


(35) 


where Tis the difference vector between the barycentres Pq and Pq. To determine 
the matrix P(si+i, Sj) it suffices to know the coordinates of four vertices in both 
coordinate systems, as well as the vector z. We can determine zhy setting P = Pq 
in eq. (35), and using the fact that T(Po) is the zero vector, because by dehnition 
the barycentre Pq is the origin of the coordinate system {x^}. Furthermore, one 
can work out by elementary trigonometry that Pq has the P-coordinates 

£'(P„) = -if'ra. (36) 


We thus obtain from eq. (35) 

= (37) 

where the last equality holds because of the identities (30).^^ Collecting all the 
information, using relations (34) and introducing the shorthand notation Xi : = 
x{Pi), T' := x'{Pi), we obtain a complete set of equations for the rotation matrix 
R{si+i,Si), namely, 

1 1 
Xi +-X4 = R{si+i,Si)xi, * = 1,2,3, X5 +-X4 = R{si+i,Si)x4. (38) 

The above construction can be generalized immediately to an arbitrary permuta¬ 
tion {** 1 , r* 2 , ** 3 , ** 4 , ** 5 } of the vertex labels {1,2, 3,4,5} for the simplex Sj, with 
{** 1 , ** 2 , ** 3 , ** 4 } denoting the vertices of the tetrahedron shared with simplex Sj+i 
and {** 1 , r* 2 , ** 3 } the vertex labels of the distinguished triangle. The relations (34), 
(36) and (37) become 


—#/ —# • 100 

- I f -I 

~ -^nii ^ ~ •^n4 ~ -^ns, -^0 ~ 2 
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The analogous relation in d dimensions is given by z = —jxiPd)- 


( 39 ) 
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Figure 3: Two-dimensional analogue of the coordinate assignments for consec¬ 
utive simplices s*, Sj+i (here represented by triangles) traversed by part of an 
oriented loop L. The vertex labels are {Pi, P2, P3} and {P[, P2, P3}, as indicated. 
The triangles share a common face a (the edge between Pi and P 2 ), and the 
common hinge around which the rotation takes place is the vertex Pi = P[. The 
coordinate system based at Pq is dehned in terms of three coordinate unit 
vectors Xi for the vertices Pi, z = 1, 2, 3, and the coordinate system {x'^} is dehned 
analogously in terms of primed quantities. The two-dimensional counterparts of 
relations (34) and (37) are given by = Xi, x '2 = T 3 and T= -x-^- 


Using a column vector notation for the coordinate four-tuples, the dehning rela¬ 
tion for the rotation matrix P can be written as the matrix equation 


(^ni T 2^^4>^’*2p 2^^4’^’35T 2^^4) Xn2) i 


(40) 

generalizing equations (38). One can check that the solution R to eq. (40) does 
not depend on the permutation of the three triangle indices and that R G 50(4). 
More specihcally, up to conjugation associated with a permutation of the hve 
vertex labels of Sj, P is equivalent to the matrix 


/I 0 0 0 \ 

0 10 0 

0 0 I ^ ’ 

Vo 0 -4^ \ j 


(41) 


which according to our discussion in Sec. 3 is an example of a simple rotation, 
in the present case representing a rotation by an angle 9 = arccos \ in the plane 
perpendicular to the triangular hinge. 

Starting at the initial simplex si, we can apply the procedure just described 
iteratively to choose vertex labels and associated coordinate systems for all other 
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simplices traversed by the loop L. However, we are not guaranteed to arrive back 
at Si with the same vertex labelling we started out with, but instead will end up 
with some label set {ai, 0 : 2 , as, a 4 , as}. The last rotation matrix in the sequence, 
R{si,Sn), therefore has to be followed by the rotation matrix Vtt implementing 
the permutation tt which brings the vertices back to the labelling chosen in si at 
the outset, that is, for which 

n{ai)=i, i = l,2,...,5, K'(^ai, ^04) = (^1, ^2, ^s, ^4) • (42) 

For the given choice of vertex labels, we have hnally arrived at the concrete 
expression for the holonomy matrix it is given by 

Rl = V^R{si,Sn) ■ ■ ■ R{s2,Si). (43) 


5.2 Choosing vertex labels (version 2) 

Any choice of coordinate systems for the simplices along a loop L will affect the 
holonomy matrix R^ at most by conjugation and the invariants (19) and (20) not 
at all. The procedure outlined in the previous subsection was natural in the sense 
that aligning the coordinate systems of successive four-simplices led to a simple 
geometric interpretation of the matrices R as being associated with the rotations 
around the triangles singled out by the loop L. However, from a computer point 
of view it is slightly inconvenient to have to identify these triangular hinges and 
to label the vertices of the four-simplices anew for each new path L, and hnally 
to compute the permutation matrix Vt^- 

It turns out to be computationally advantageous to dehne local coordinate 
systems which make use of the fact that as part of the Monte Carlo set-up each 
four-simplex already comes with a labelling of its vertices in terms of a permuta¬ 
tion of the indices {1,2,3,4, 5}. We have to generalize our considerations of the 
previous section only slightly to obtain the holonomy matrix for this case. 

Consider again the transition from simplex Si to simplex Sj+i. The general 
situation is that the four vertices spanning the common tetrahedral face between 
the two four-simplices have labels (ui, n 2 , ns, 77 , 4 } in s, and labels {ki, k 2 , k^, k^} 
in Sj+i. Let and k^ be the remaining labels of the hfth vertex in Sj and Sj+i 
respectively. The transformation matrix i?(sj+i,Sj) we are looking for satishes 


(^fcl T 2^^5 ’ 2^^5’^7i;4T 2^^^ .^('^ 7 + 1 7 ^^ 7 ) ‘ (^ni 5 ^n 2 5 ^713 7 ^n 4 ) 5 

(44) 

in analogy with eq. (40). From this equation we can read off that the solution 

Ri^^i+l j T 2 ^^5 ’ ^^2 T 2 ^^5 ’ ^^3 4 " 2 ^^5 ’ ^^4 2 ' (^fii) ^n 2 ) j 

(45) 
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is invariant under simultaneous permutations of the label sets {ni, ^ 2 , ns, n 4 } and 

{ki,k2,h,ki}. 

Note that in the computer program no attention is paid to the relative ori¬ 
entation of the four-simplices when labelling their vertices with a permutation 
of {1,2, 3,4, 5}. This implies that a matrix computed from (45) can 

have determinant —1 rather than 1, i.e. belong to 0(4) rather than 5*0(4), which 
according to Sec. 3 is the generic case anyway. However, for a closed curve L 
there will always be an even number of 0(4)-matrices with determinant —1 and 
the holonomy matrix Rl will therefore always be a S'0(4)-matrix. In terms of 
the new /^-matrices determined from (44) we can write the holonomy as 

R'^ = R{si,Sn)---R{s2,si). (46) 

This is the analogue of (43), but without the need for a permutation matrix 
since we will always return to the same labelling of the vertices of the simplex 
Si where the loop L starts and ends. We will use the prescription leading to 
expression (46) in the measurement of holonomies described in Sec. 7 below. 

By construction only a hnite number of i7-matrices can occur. They can be 
computed, stored and looked up for given pairs of vertex label sets, without having 
to perform the matrix inversion and multiplication of formula (45) each time. In 
addition, we can make maximal use of the permutation invariance mentioned 
above. For example, we only need the hve matrices 



/ tti /3i 7i l\ 

q;i l3i 471 0 

tti 3/^1 0 0 

^^ 2 q;i 0 0 Oy 


/ -ai /3i 7i l\ 

—«! /3i 471 0 

—«! 3/^1 0 0 ’ 

\- 2 ai 0 0 0 / 


(47) 


/ 0 -2/3i 7 i 1 \ / 0 0 -371 1 \ / 0 0 371 - 1 \ 

0 -2/3i 471 0 0 2/3i -471 0 0 2 / 3 i -71 -1 

—3/^1 0 0’ tti —/3i —471 0 ’ tti —/3i —71 —1 

\-tti -3/3i 0 0 / \-ai -I3i -471 0 / \-ai -71 - 1 / 

with 



when determining i?(sj+i, Sj), because we can always hud a simultaneous permu¬ 
tation of (ni, n 2 , ns, n 4 } and {ki, k 2 , k^, k^} to make the second matrix on the 
right-hand side of eq. (45) match one of these hve matrices. 


6 Wilson lines in CDT 

Having given a concrete prescription for the computation of the transition matri¬ 
ces i?(sj+i, Sj) and how to obtain from them the holonomy Rl associated with a 
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given lattice loop L, let us now turn to what we can learn about the behaviour 
of the holonomies and their potential use in probing “quantum geometry”, that 
is, the geometric properties of the dynamically generated ground state of the 
nonperturbative path integral. 

One possibility would be to study the holonomy group of a given, hxed piece- 
wise flat geometry, taken to be a typical member of the gravitational path integral 
ensemble generated by the Monte Carlo simulations. We know that these geome¬ 
tries are non-differentiable and highly singular (similar to the conhgurations - 
the “paths” - of an ordinary non-relativistic quantum-mechanical path integral), 
and one could compare the holonomy groups of different, fixed such geometries 
with those of smooth Riemannian manifolds. 

We will be interested here in the physically more interesting case of us¬ 
ing holonomies to construct observables in the fully dynamical, nonperturbative 
quantum theory. By “observables” we mean in this context coordinate-invariant 
quantities, which are operationally well-dehned on the quantum-fluctuating en¬ 
semble of geometries. We do not require them to be related explicitly to any 
truly observable phenomenological effects (other than perhaps in some semiclas- 
sical limit), which would be a tall order in any theory of quantum gravity. An 
example of an observable in this looser sense is the (expectation value of the) 
spectral dimension of quantum spacetime, a quantity which has been measured 
explicitly in CDT quantum gravity [8], and also studied in other formulations 
[19]. 

Coming up with physically interesting observables in the sense just described 
is still a formidable challenge in background-free, nonperturbative quantum grav¬ 
ity. To illustrate the point, consider some two-point function G 2 {x,y) in standard 
quantum held theory on a fixed background. Its naive analogue on a nonpertur- 
batively huctuating ensemble of geometries is not a well-dehned observable, be¬ 
cause there is no coordinate-invariant way to refer to the same two points x and 
y throughout the ensemble. One workaround is to specify the geodesic distance 
of the two points to be r and integrate over all possible positions of x and y 
subject to this constraint. Averaging this quantity in the path integral over the 
geometric ensemble then leads to a well-dehned two-point function 

Q2{r)=jV[gf,^]e-^^^^''^Jdx dy^g{x)g{y)G2{x,y)6{r-dg^Xx,y)) (49) 

(see [20] for a concrete implementation in two dimensions). The dependence of 
the diheomorphism-invariant propagator Q 2 on the coordinate-invariant geodesic 
distance r captures nontrivial physical information. 

There are similar issues when trying to construct an observable that depends 
not just on two points, but on an entire closed curve in spacetime. An obvious 
generalization of the prescription leading to a well-dehned two-point function 
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would be to consider the trace invariants (19), (20) for a subclass of loops sharing 
certain invariant geometric features in terms of their length and shape, and then 
to integrate over all possible locations of such loops. 

In the present work, we will pursue a similar strategy to construct and measure 
a particular class of well-dehned Wilson loop observables, but instead of referring 
to intrinsic geometric properties of the underlying paths, we will introduce a 
dynamical point particle, couple it to the quantum geometry, and compute the 
holonomy along its world line. In other words, we will consider an interacting 
system of matter and geometry, where a massive point particle is coupled to 
pure quantum gravity, given in terms of the usual CDT ensemble of fluctuating 
geometries. The (Euclidean) action of a point particle of mass m with spacetime 
trajectory 7 in a Riemannian geometry with metric is given by the mass times 
the proper length of 7 , 

^P-P- = mjdl, dl= ^ds. (50) 

Concretely, we will add an appropriate simplicial version of the point particle 
action (50) to the Einstein-Hilbert action and update the combined, interacting 
system using Monte Carlo simulations, as described in Sec. 7 below. 

As usual in CDT, we consider spacetimes of topology x , where - for 
convenience - time has been compactihed to a circle. The particle world lines 
whose holonomies we will measure are oriented in positive time direction and 
wind around the periodic time direction exactly once. Accordingly, our loops all 
have macroscopic length, with a minimum that depends on the time extension 
of the spacetime. In line with standard terminology from gauge theory we will 
refer to these holonomies and their associated trace invariants (the context should 
make clear which is meant) as Wilson lines. 

To understand the geometry of these Wilson lines, we need to recap briefly 
some aspects of the triangulated spacetimes in the conhguration space of the 
path integral. More complete descriptions can be found elsewhere, see [5] and 
references therein. In standard CDT, each spacetime T has a proper-time slicing 
with integer label t, and is assembled from four-simplices in a layered fashion^^, 
where one layer of thickness At = 1 is a piecewise flat piece of spacetime of 
topology X /, all of whose vertices are contained in either of its spatial boundary 
submanifolds at times t or t-|- 1. These submanifolds are arbitrary triangulations 
in terms of equilateral tetrahedra, and all have the topology of a three-sphere. An 
entire four-geometry of proper-time extension T is obtained by gluing together 
T subsequent layers along matching three-geometries, and hnally identifying the 
hnal boundary of the last layer with the initial boundary of the hrst layer. 

^^see [21] for a generalization of CDT geometries, without strict time slicing, but maintaining 
causality 
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(4,1) (3,2) 

Figure 4: The two types of elementary simplicial building blocks used in CDT 
quantum gravity, and how they are positioned with respect to two adjacent slices 
of constant proper time. Together with the (4,1)- and the (3,2)-simplex shown, 
also their time-reversed versions, the (1,4)- and the (2,3)-simplex, appear. 


Before the Wick rotation, we distinguish between space- and timelike edges. 
The former are always contained in a three-dimensional spatial slice of hxed 
proper time, whereas the latter interpolate between two adjacent spatial slices 
of hxed proper times t and t + 1. There are (up to time rehection) two types 
of elementary building blocks, the (4,l)-simplex and the (3,2)-simplex, see Fig. 
4. The notation indicates how they are positioned in a given layer with respect 
to the slices of constant integer time: a four-simplex of type (i, j) has i vertices 
with time label t and j vertices with time label t + 1. Building blocks of different 
types have different numbers of time- and spacelike edges. Although after the 
Wick rotation all edges are spacelike and - with the convention adopted in this 
paper - of equal length £, one can still distinguish the different building blocks 
because of how they are positioned inside the layered structure, since the Wick 
rotation leaves the connectivity of the triangulation intact. 

7 Measurements 

As usual in CDT simulations, we keep the number T of proper-time steps hxed, as 
well as the total four-volume (dehned as the total number of four-simplices), 
which turns out to be technically convenient. For the measurements reported 
below, we worked with T = 80 and = 20.000, at the point (kq. A) = (2.2, 0.6) 
in the space of coupling constants, which lies inside the “de Sitter phase” of CDT 
quantum gravity [3] (kq is the bare inverse Newton constant, and A denotes the 
so-called asymmetry parameter). The action we use is the Regge form of the 
Einstein-Hilbert action [17], applicable to piecewise linear geometries. On causal 
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dynamical triangulations with identical equilateral building blocks it assumes a 
particularly simple form (see [5] for details). 

The Monte Carlo computer simulations proceed from an initial conhguration 
put in by hand, after which we apply a standard Metropolis algorithm to update 
the geometry, using as action the Regge-Einstein-Hilbert action, and allowing 
only those Monte Carlo moves (changes in the geometry) that are compatible with 
the foliated structure of the spacetimes. The simulations generate a sequence of 
four-dimensional geometries - spacetime histories represented by triangulations - 
which after sufficiently many updates will be independent of the chosen starting 
conhguration. Computing the expectation value of a given observable amounts 
to measuring the observable many times for statistically independent geometries 
generated by the Monte Carlo simulation, and calculating its average inside the 
path integral. 

We describe next how to add a Wilson line to the Monte Carlo simulations of 
pure gravity. In the selected initial conhguration, we put in by hand a closed path 
that winds once around the periodic time direction. Without loss of generality, 
we follow the choice made in Sec. 3 and consider only piecewise straight paths 
connecting the centres of neighbouring four-simplices. It is important to under¬ 
stand that not every type of four-simplex can be a neighbour of every other type 
of four-simplex, in the sense of having a three-simplex in common. The reason 
for this is that before the Wick rotation the three-simplices also come in diherent 
types, depending on the time- and spacelike character of their edges, and that 
gluings of two four-simplices along two tetrahedral faces are only possible if the 
metric properties of the tetrahedra match exactly. For this reason, a (4,l)-simplex 
cannot be a neighbour of a (2,3)-simplex, say. 

For a path moving forward in time this implies that only particular sequences 
of the simplex types can occur along it, namely, those of the form 


. (4,1),..., (4,1), (3, 2),..., (3,2), (2, 3),..., (2, 3), (1,4),. 

^-V-^ ^-V-^ ^-V-^ ^- 


/I A\ (A 


(51) 

which should be continued cyclically and read as: a (4,l)-simplex can only be 
followed by another (4,l)-simplex or by a (3,2)-simplex, a (3,2)-simplex can only 
be followed by another (3,2)-simplex or by a (2,3)-simplex, and so forth. It follows 
that such a path needs to go through at least four four-simplices to pass from 
one layer to the next, and to arrive at the same type of four-simplex it started 
from. For example, consider a path starting at the barycentre of a (4,l)-simplex 
in the layer between times t and t + 1. It has to pass through at least three other 
four-simplices in the same layer before arriving at a (4,l)-simplex in the layer 
between times t + 1 and t + 2, namely, one (3,2)-simplex, one (2,3)-simplex and 
one (l,4)-simplex. In other words, at least four steps are necessary to advance by 
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one time unit At = 

In our set-up, an oriented path associated with a Wilson line is not allowed 
to go back in time relative to the proper-time foliation and the sequence dehned 
in (51), thereby enforcing some degree of “causality”. For example, it may move 
among the (3,2)-simphces of a given layer, but not subsequently go back to a (4,1)- 
simplex of the same layer; it can only proceed to a (2,3)-simplex, as specihed by 
(51). Another restriction put on the paths is that they are not allowed to self¬ 
intersect, mainly to prevent them from meandering for very long times inside a 
given layer. 

The Boltzmann weights of the combined configurations of four-geometry and 
particle path now contain an additive contribution to the Einstein-Hilbert action, 
which is the lattice version of the continuum action (50) and simply given by 

A/-P- = moNL, (52) 

where tuq is the bare particle mass and the number of four-simplices en¬ 
countered by the loop L. During the Monte Carlo simulations, the path evolves 
in computer time along with the geometry. The usual local update moves are 
performed on the geometry. Whenever the particle world line happens to pass 
through simplices affected by such a local rearrangement, it will be broken up 
there. Part of the updating algorithm is then to determine all possible ways in 
which the two loose ends of the path can be reconnected, and weighing them with 
the appropriate Boltzmann weights. Lastly, moves will be disallowed if they lead 
to a path which self-intersects or goes backward in time.^^ 

The behaviour of the particle path L depends on the parameter rriQ in the 
action For large mo, long paths are strongly suppressed, and the length will 

be Nl = 4T = 360, which minimizes the action (52). The constraints forbidding 
backtracking in time and self-intersections of L are in this case irrelevant, because 
such conhgurations occur only very rarely in the computer simulations. As mo 
becomes smaller, longer loops become less costly and the length of L will become 
longer and fluctuate more. As a consequence, the constraints play a nontrivial 
role in limiting the length of the loop. 

Once the simulations are running and have thermalized, we sample conhgura¬ 
tions, measure Rl along the particle world line and compute the trace invariants 
ti and t 2 from (19) and (20). This allows us to determine the angles 6 i and 

^^Because of this substructure, one can make a further subdivision of time, with units of 1/4, 
which is sometimes useful. This was first introduced in [7], see also [5]. 

^"‘We will describe the Monte Carlo simulations of a point particle coupled to CDT quantum 
gravity in detail elsewhere, where we also investigate systematically the effect of the particle on 
the quantum geometry and vice versa. In the present article we focus on the construction and 
measurement of the Wilson lines associated with the particle paths, keeping the account of the 
technicalities of the computer simulations to a minimum. 
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Figure 5: Histogram of the invariant angles {61,62) characterizing a SO(4)- 
holonomy Rl, from Monte Carlo measurements. 


62 in the interval [0, 2tt\ up to an interchange 61 o 62 and up to the reflections 
6i —)■ 2 n — 6i. Fig. 5 shows the histogram of the measured values for the angles 
6i. Its shape is independent of the mass mo in the point particle action (52), and 
by construction displays the above exchange and reflection symmetries. We have 
verihed that after normalization the measured distribution is in perfect agree¬ 
ment with the theoretical distribution P{6i, 62) of (26), which was derived under 
the assumption that the holonomy Rl is uniformly distributed over the group 
manifold ^(^(d). 

8 Discussion and outlook 

One could have wondered a priori whether the coordinate-free set-up of (Causal) 
Dynamical Triangulations is suited to describing holonomies and Wilson loops. In 
this article, we have demonstrated conclusively that it is straightforward to dehne 
and compute these quantities. More than that, we hnd it difficult to envisage 
a framework, lattice-based or otherwise, that would make the computation of 
Wilson loops on four-dimensional curved manifolds even simpler. Our explicit 
construction involved a particular choice of Cartesian coordinate systems on the 
individual simplicial building blocks, but clearly many other choices are possible 
and would not affect the hnal result, which was formulated in terms of coordinate- 
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independent quantities. 

One could also have worried a priori that the computation of holonomies in 
the (C)DT framework was affected strongly by discretization effects, especially 
since the equilateral simplices have just a single interior angle ctint = arccosl /4 
(the angle between two three-dimensional faces sharing a two-dimensional hinge), 
despite the fact that this angle is irrational. In addition, we saw in Sec. 5 that only 
a hnite number of different transition matrices R{si^i, Si) occur in the holonomy 
computations. However, we have not observed any sign of such discretization 
artefacts. On the contrary, our main result is that for the class of Wilson lines 
considered, the holonomies appear to cover SO( 4 ) densely and uniformly. This 
also implies that the holonomy group of the quantum geometry generated in the 
de Sitter phase of CDT quantum gravity is SO( 4 ), which is the same as the 
holonomy group of a generic orientable four-dimensional Riemannian manifold. 
Conversely, we have not found any tendency of our macroscopic Wilson lines to 
cluster around the identity of the group SO( 4 ). 

We have shown that the CDT framework is well suited for investigating Wil¬ 
son loops in nonperturbative quantum gravity. Our Wilson lines are well-dehned 
observables, but we have at this stage no direct physical interpretation to re¬ 
late them to specihc classical or quantum properties of the underlying quantum 
spacetime. The challenge for any theory of quantum gravity is to come up with 
observables which do this. As already mentioned in the Introduction, quantities 
involving holonomies in one way or other are natural candidates for encoding 
information about the curvature of (quantum) spacetime. Our analysis of the 
holonomy of minimal loops around a single triangular hinge shows that the clas¬ 
sical relation between the holonomy of an inhnitesimal planar loop and the local 
curvature - here in the form of a dehcit angle a la Regge - continues to hold on 
the piecewise flat geometries of (Causal) Dynamical Triangulations. 

From the point of view of the regularized lattice formulation, these minimal 
loops are not particularly interesting, since they merely probe geometry at the 
cutoff scale, which is dominated by lattice artefacts, that is, the largely arbi¬ 
trary details of the regularized set-up at this scale. On the other hand, the large 
gravitational Wilson loops we have studied do not obviously contain retrievable 
curvature or other geometric information, unless it is hidden in higher-order corre¬ 
lators. Here one should of course keep in mind that even in the classical continuum 
theory we do not know how to relate the values of non-inhnitesimal Wilson loops 
on a general Riemannian manifold to its curvature, because of the nonabelian 
nature of the metric connection. The most obvious quantities to try to dehne 
and investigate are therefore lattice Wilson loops which are much bigger than 
the minimal loops, but sufficiently small to have an interpretation in terms of a 
suitably averaged curvature, in the continuum limit. The challenge is at least 
two-fold: (i) to dehne suitable classes of closed curves which have an invariant 
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meaning when we integrate over all geometries, and (ii) to use them to hnd a 
notion of (quantum) curvature whose expectation value remains hnite and well- 
dehned at the Planck scale, while converging to some function of the continuum 
Riemann tensor in the classical limit. The Regge dehnition of (scalar) curvature 
on piecewise linear geometries is simple, but offers little insight into how to dehne 
a well-behaved notion of quantum curvature. The use of holonomies may offer an 
alternative, more flexible tool to achieving this goal. Further research into this 
issue is currently under way. 

Acknowledgments 

JA and AG acknowledge support from the ERG Advanced Grant 291092 “Ex¬ 
ploring the Quantum Universe” (EQU) and by FNU, the Free Danish Research 
Gouncil, through the grant “Quantum Gravity and the Role of Black Holes”. 
JJ acknowledges the support of grant DEG-2012/06/A/ST2/00389 from the Na¬ 
tional Science Gentre Poland. The contribution of RL is part of the research 
programme of the Foundation for Fundamental Research on Matter (FOM), hnan- 
cially supported by the Netherlands Organisation for Scientific Research (NWO). 


References 

[1] J. Ambjprn, J. Jurkiewicz and R. Loll, Phys. Rev. Lett. 85 (2000) 924-927 
[hep-th^002050]; Nucl. Phys. B 610 (2001) 347-382 [hep-th/0105267]. 

[2] J. Ambjprn, J. Jurkiewicz and R. Loll, Phys. Rev. Lett. 93 (2004) 131301 
[hep-th^404156]; Phys. Rev. D 72 (2005) 064014 [hep-th/0505154]. 

[3] J. Ambjprn, A. Gorlich, J. Jurkiewicz and R. Loll, Phys. Rev. Lett. 100 
(2008) 091304, arXiv:0712.2485 [hep-th]; Phys. Rev. D 78 (2008) 063544, 
arXiv:0807.4481 [hep-th]. 

[4] J. Ambjprn, J. Jurkiewicz and R. Loll, in Approaches to quantum gravity, 
ed. D. Oriti (Gambridge University Press, Gambridge, UK, 2009) 341-359 
[hep-th/0604212]; Annalen Phys. 19 (2010) 186-195; in Foundations of space 
and time, eds. G. Ellis, J. Murugan and A. Weltman (Gambridge University 
Press, Gambridge, UK, 2012) 321-337, arXiv: 1004.0352 [hep-th]; PoS LAT- 
TIGE 2010 (2010) 014, arXiv:1105.5582 [hep-lat]. R. Loll, Glass. Quant. 
Grav. 25 (2008) 114006, arXiv:0711.0273 [gr-qc]; J. Ambjprn, A. Gorlich, 
J. Jurkiewicz and R. Loll, arXiv:1007.2560 [hep-th]; in Ftandbook of Space- 
time, eds. A. Ashtekar and V.Petkov (Springer, Heidelberg, 2014) 723- 
741, arXiv:1302.2173 [hep-th]; Int. J. Mod. Phys. D 22 (2013) 1330019, 



arXiv: 1305.6680 [gr-qc]; J. Ambjqrn, S. Jordan, J. Jurkiewicz and R. Loll, 
AIP Conf. Proc. 1514 (2012) 60-66, arXiv:1302.2181 [liep-th], 

[5] J. Ambjqrn, A. Gorlicli, J. Jurkiewicz and R. Loll, Phys. Rept. 519 (2012) 
127-210, arXiv: 1203.3591 [liep-th]. 

[6] J. Ambjprn, A. Gbrlich, J. Jurkiewicz, A. Kreienbuehl and R. Loll, Glass. 
Quant. Grav. 31 (2014) 165003, arXiv:1405.4585 [hep-th]. 

[7] J. Ambjprn, A. Gorlicli, J. Jurkiewicz, R. Loll, J. Gizbert-Studnicki and 
T. Trzesniewski, Nucl. Phys. B 849 (2011) 144-165, arXiv:1102.3929 [hep- 
th]. 

[8] J. Ambjprn, J. Jurkiewicz and R. Loll, Phys. Rev. Lett. 95 (2005) 171301 
[hep-th/0505113]. 

[9] R. Williams, in Approaches to quantum gravity, ed. D. Oriti (Gambridge 
University Press, Gambridge, UK, 2009) 360-377. 

[10] J. Ambjprn and J. Jurkiewicz, Phys. Lett. B 278 (1992) 42-50; M.E. Ag- 
ishtein and A.A. Migdal, Mod. Phys. Lett. A 7 (1992) 1039-1062. 

[11] K.G. Wilson, Phys. Rev. D 10 (1974) 2445-2459. 

[12] G. Modanese, Phys. Rev. D 49 (1994) 6534-6542 [hep-th/9307148]. 

[13] D. Diakonov and V. Petrov, J. Exp. Theor. Phys. 92 (2001) 905-920 [hep- 
th/0008035]. 

[14] T. Thiemann, in Approaches to quantum gravity, ed. D. Oriti (Gambridge 
University Press, Gambridge, UK, 2009) 236-252; G. Rovelli, Glass. Quant. 
Grav. 28 (2011) 153002, arXiv:1012.4707 [gr-qc]. 

[15] G. Glarkson, G. Ellis, J. Larena and O. Umeh, Rept. Prog. Phys. 74 (2011) 
112901, arXiv:1109.2314 [astro-ph.GO]. 

[16] H.W. Hamber and R.M. Williams, Phys. Rev. D 76 (2007) 084008, 
arXiv:0706.2342 [hep-th]; Phys. Rev. D 81 (2010) 084048, arXiv:0907.2652 
[hep-th]. 

[17] T. Regge, Nuovo Gim. 19 (1961) 558-571. 

[18] B. Simon: Representations of hnite and compact groups. Graduate Studies 
in Mathematics, vol. 10 (American Mathematical Society, 1996). 


29 



[19] O. Lauscher and M. Reuter, JHEP 0510 (2005) 050 [hep-th/0508202]; 
P. Hofava, Phys. Rev. Lett. 102 (2009) 161301, arXiv:0902.3657 [hep- 
thj; M. Reuter and F. Saueressig, New J. Phys. 14 (2012) 055022, 
arXiv: 1202.2274 [hep-thj. 

[20] J. Ambiprn and K.N. Anagnostopoulos, Nucl. Phys. B 497 (1997) 445-478 
[hep-lat/9701006]. 

[21] S. Jordan and R. Loll, Phys. Lett. B 724 (2013) 155-159, arXiv:1305.4582 
[hep-th]; Phys. Rev. D 88 (2013) 044055, arXiv: 1307.5469 [hep-th]. 


30 



